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ABSTRACT 

We propose a random matrix model as a representation for D = 1 open strings. 
We show that the model is equivalent to iV fermions with spin in a central potential 
that also includes a long-range ferromagnetic interaction between the fermions that 
falls off as l/(rjj) 2 . We find two interesting scaling limits and calculate the free 
energy for both situations. One limit corresponds to Dirichlet boundary conditions 
for the dual graphs and the other corresponds to Neumann conditions. We compute 
the boundary cosmological constant and show that it is of order 1/ log(/3). We also 
briefly discuss a possible analog of the Das-Jevicki field for the open string tachyon. 
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Matrix models have lead to many insights in low dimensional string theory. 
The model that is perhaps the simplest to understand, yet ironically has the richest 
structure, is that which corresponds to a scalar field theory on a line coupled to 
gravity. This is the matrix quantum mechanics model [1] with the matrix potential 
tuned such that the theory is almost critical [2,3-5]. This model is equivalent to a 
theory of closed strings in 1 + 1 dimensions, where time is one coordinate and the 
eigenvalues of the matrices A, are related to the Liouville coordinate [6-9]. 

A natural extension of this work is to find the corresponding theory of open 
and closed strings. One suggestion is to add the term £Trlog(0(t) — /i) to the 
matrix potential, and tune £ and [i to some appropriate values [10]. The trouble 
with this is that, unlike the zero dimensional case [11-13], this term does not arise 
by integrating out fields which are 1 x N matrices that couple to <p. These fields are 
fundamental representations of a global U (N) group and it is precisely these fields 
that generate boundary terms in the Feynman diagrams. The problem with the 
log term is that the kinetic piece of the fundamental fields has been ignored. One 
might argue that it is all right to drop this term if the mass and the couplings are 
very large. However, it has been pointed out that the interesting critical behavior 
occurs when the argument of the log approaches zero [10,14]. This is precisely 
where the kinetic piece should become important. 

In this paper we explicitly include these fundamental fields in the full la- 
grangian. We demonstrate that the theory is equivalent to N fermions with spin 
which interact among each other with long-range ferromagnetic interactions. We 
show that there are two interesting scaling limits for the masses of the fundamen- 
tal fields. One limit has the mass diverging logarithmically as iV — > oo, while the 
other case has a mass that diverges, but at a milder rate. The former case appears 
to correspond to an open string theory with Dirichlet boundary conditions on the 
dual graphs (and hence Neumann conditions on the original graphs), while the 
latter limit corresponds to a theory with Neumann conditions. In both cases we 
compute the free energy as an expansion in the open and closed string couplings. 
The Neumann case has almost all terms in the expansion negative definite, while 
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the Dirichlet case has the same form for the free energy, but with the opposite sign 
for the open string coupling. The Neumann case also has a critical value for the 
open string coupling, at which point divergences start appearing in the free energy. 
We also show that the boundary cosmological constant scales as (log/9) -1 . Finally, 
we discuss a possible analog to the Das-Jevicki collective coordinate field [6] for 
the open string tachyon. 

To begin, consider a triangulated surface £, with boundary B, where B is not 
necessarily connected. Suppose that U is the world-sheet for a string theory which 
has a coordinate t. Each vertex of a triangle marks a particular point in t space, 
t v . Then an edge, £, connecting two vertices has associated with it a difference in 
t, Ate — t Vl — t V2 . The complete partition function should be comprised of a sum 
over all assignments of t v to the vertices. This is equivalent to summing over all 
Ate, so long as we impose the constraint 



where the sum over £ is over the three edges of any triangle C. 

The string world sheets should also have boundary conditions. For the trian- 
gulated surface S, Neumann boundary conditions impose the constraint Ate = if 
£ is an edge that intersects the boundary. On the other hand, Dirichlet boundary 
conditions impose the same constraint Ate — 0, but now £ is an edge that lives on 
the boundary. 

The dual surface of E is found by bisecting all edges. This operation maps 
vertices to faces and vice versa, while edges are mapped to new edges at 90 degrees. 
The dual surface can be generated by a field theory of iV x N matrices, a ^, with 
cubic interactions [15,2]. The surfaces are the Feynman diagrams for this theory, 
with the propagators forming the edges. Ate f° r an edge in E maps to pi, the 
momentum flowing through the bisecting propagator, while the constraint in (1) 
naturally maps to the constraint that the sum of the momenta entering any vertex 
is zero. 




(1) 
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The boundaries of E present a slight problem when finding the dual graph. 
The edge dual to a boundary edge has nothing to attach to outside the boundary. 
We can remedy this by introducing new fields, ip a and Xa which transform in 
the fundamental representation of a global U(N) and which couple to a ^. The 
propagators of these fields can then be used to tie off the ends of the propagators 
dangling over the boundary. ip a and Xa ar e assumed to be fermionic and we will 
see later that two fields are necessary in order to have a nontrivial theory. 

The boundary conditions in U lead to analogous constraints for the dual 
graphs [16]. Edges that are normal to the boundary are mapped to edges that 
are parallel. Therefore, Neumann conditions on U correspond to dual graphs with 
zero momentum in the (f) propagators that run along the boundary. This means 
that if a large amount of momentum flows into a region of the Feynman diagram 
which includes a long stretch of the boundary, then the momentum should mainly 
flow through the propagators of the fundamental field and not the propagators 
that are next to but not actually intersecting the ip or x propagators. This is 
accomplished by making the masses of ip a and Xa very large, in which case, the 
weights of the Feynman diagrams strongly favor a large amount of momentum 
flowing through the ip and x propagators. As the masses diverge, the correlation 
lengths of these fields shrink to zero and all fermion loops become localized in t. 
Thus, Neumann conditions on U lead to Dirichlet conditions for the dual theory. 
Likewise, Dirichlet conditions on £ lead to dual graphs with zero momenturm 
flowing off the boundaries. Such graphs will dominate the free energy if the masses 
of ipa and Xa are chosen to be small. The dual graph then has Neumann boundary 
conditions. 

Let us consider the action, previously considered by Yang [14] 



We have chosen to couple the fermion fields to the square of the hermitian matrix 



S = /3 dt\ Tr ( -0 2 - V(<f>)) + i^Ja + ixlt 




a 



(2) 
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for later convenience, but this won't affect the critical behavior in the double scaling 
limit. It is necessary to couple x to the complex conjugate of in order to build 
open string states that are invariant under a global SU (N) rotation. 

In the large N limit the free energy derived from this action is dominated by 
the planar diagrams and can be expressed as an expansion in 1/N, 

F = - N 2 ^N- h (N/f3) &xe& {-i 2 N/f3) len & h F{Z,g,h^), (3) 

surfaces 

where g is the genus of the surface and h is the number of holes. Each vertex 
in the dual diagram is weighted by (N/(3) l l 2 and each vertex on the boundary 
is weighted by ^(N/fl) 1 / 2 . Hence in terms of the original triangulated graph, 
these weights lead to the terms in the free energy exponentiated by the area and 
length in appropriate units. J-{£, g,h, /i) is a weight for the Feynman diagram 
corresponding to a particular surface. If \i is very large then the \i dependence of 
T is approximately (l//i) length . 

To proceed, let us diagonalize the matrix 0, = WAU, where A is diagonal and 
U is an element of SU(N). Moreover, let us rotate the fermion fields ip a — > u\ b i>b 
and Xa — > XbUba- Then the action (2) becomes 

S = f3 J dtj^(h 2 a - V(X a ) + i^j a + IxlXa 

+ li'l^a + IxWaXa - ^\^a ~ HxlXa) 

- E iypu\ h {\ a - x b ) 2 (uu^) ba + tyKurtUik - i{UU\ bX \xa) ■ 

(4) 

Letting A = iUU\ we see that the action is quadratic in A and hence can be inte- 
grated out after shifting variables [17,18]. To this end, consider the path integral 
for propagation of the vacuum state from time t — to time t — T, 



Z = J V<\)Vi)V^Vx^Vxe iS 

= I V\VUV^V4>Vx ] VxY[(A(\(t)) 2 e iS , 

J 4- 



(5) 
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where A(A(t)) is the vandermonde determinant 

A(X(t)) = H(\ a (t)-X b (t)f. (6) 

a<b 

Completing the square and shifting A, leads to the extra term in the action [18] 

J * 5 — — • ( ' 

Noting that U\t + At)U(t) = 1 + iA(t + At/ 2) At + we can replace the U 
integrals with A integrals, less one U integration because of the global SU(N) 
symmetry. Since the diagonal elements of SU(N) commute with A, we should 
restrict the integration to a subspace of A. If we only consider the insertions of 
SU(N) invariant states into the path integral, then this subspace is just the off- 
diagonal elements of A. This integration over A ab leads to factors that cancel off 
the vandermonde determinants in the Jacobian, except for a leftover contribution 
at each end point because of the extra integration over an SU (N) variable. Hence 
the partition function reduces to 

Z = J PA^tp^ x tp x A(A(0))A(A(T))e iScff , (8) 

where S e s is the effective action after integrating out A. 

The fermion kinetic terms lead to the standard quantization for ip a and Xa, 

{>a, = Sab, {Xa, x\} = Sab, (9) 

where a factor of /9 1//2 has been absorbed into the fermion fields. Therefore, the 
hamiltonian for SU (N) invariant states is given by 

H = E (~2^w a + v{Xa) + ^° A M + *° A 5x£) - + XaX ^) 

t t t t ( 10 ) 

1 \- (Wb - X b Xa)(%^ ~ *M 

The last term in the hamiltonian leads to a repulsive force between the eigenvalues. 
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However, because ip a and Xa don't commute with tjr a and xl, respectively, it 
is necessary to determine the ordering of the operators in H so that the theory is 
well defined. The ordering follows from the definition of the time derivative in the 
path integral. As is the case with the kinetic term for the hermitian matrix, the 
kinetic term for the fermions should link an otherwise uncoupled chain of fields. 
With this in mind, let us define the kinetic term for ip a as 

t 

Rotating ip a by and integrating out the angular variables leaves the following 
contribution to the path integral (ignoring the x field): 

z^x) = / v^v^^ v UY}^^l{t){Mt) - W - At)) 

J ^ t=o a 

+ At(^ut)\ 2 a (t)Mt) - ^im a (t))] (12) 

v At 4>l(t)ipl(t)Mt - &t)Mt - At) 

(A a - \ b ) 2 

where the integrations over the variables are at the points between t = and 
t = T inclusive. Then, with no tp insertions in Z^(X), (12) reduces to 

Z^X) =Cn(l-^At^(A2(t)-/i)) 

t a 
T 

=Cexp(-t/? J <ft^(*)-4 



(13) 



where C is an unimportant constant. If we now insert the SU(N) invariant operator 

c aia2...ajv 

ipl 1 ipl 2 ...ipl N into Z^x) at t = 0, we find that the contribution to the path 
integral is C. This means that there is no contribution to the energy from the 
four fermi interaction if the ip states are either completely full or empty. Thus we 
learn that the explicit ordering given by the hamiltonian in (10) follows from our 
definition of the infinite matrix chain. 
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From the hamiltonian (10) we can determine the ground state, which should 
be a state that is SU(N) invariant. Given an SU(N) invariant state other 
invariant states can be constructed by acting on |\&) with not only the operator 
tr0 n , but also with J™, J! 1 and Jl 1 , where these last three operators are defined as 

The J n operators form an algebra 

\J n z , J%] = ±Jl +m [Jl, J m ] = 2J? +m , (15) 

and in particular, the operators J = are the generators of an SU (2) lie alge- 
bra. Furthermore, after diagonalizing <p and redefining ip and \ we can construct 
operators S+ a , S- a and S za , where 

S+a=lf>lXa> S- a = Xa^a S za = ^(^ a ~ Xaxl) , (16) 

— * — * — * 

and which satisfy the relation J n = J2 a ^a^a- Clearly, for each index a, S a can be 
interpreted as a spin operator. Hence the hamiltonian (10) can be reexpressed as 

"=E(-^i + v- ( Aj-|(^-i)(^-,)) + ^E ! f^- (») 

The result is a system equivalent to N fermions with spin in a central potential 
which includes a position dependent magnetic field and with Heisenberg ferromag- 

— * — * 

netic long range couplings between the spins. Although J ■ J does not commute 
with H, J z does, so the energy eigenstates can be classified by their J z quantum 
numbers. 

It is not known how to calculate the exact spectrum of (17). However, for 
some choices of 7 and ji we can at least find the ground state. First consider the 
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case 7 < 0. If [i > 0, then the magnetic field, while position dependent, points 
in the same direction for all values of A a . Hence, the ground state | j), nas a ^ 
spins pointing down in the direction of the field and J z = —N/2. The reduced 
hamiltonian for spins with this configuration is [14] 

*i=e(-^J+^>+|(^-4 (is) 

Since all spins are down, the fermion interaction term has dropped out. If we now 
let /j, < 0, then | j) is not necessarily the ground state anymore, since the field 
is no longer pointing in the same direction for all A. For some negative value of 
/i, the lowest energy state with J z = —N/2 becomes degenerate with the lowest 
energy state with J z = 1 — N/2. Hence, this value of /i is a critical point for the 
fermion mass. 

Next consider the case 7 > 0. If yu < 0, then the field is pointing in the same 
direction for all A, but now the ground state has all spins up. Examining the 
hamiltonian (17), we find that there is no longer any 7 dependence in the reduced 
hamiltonian. However, if we assume that V(X) is infinite in the range |A| > a, 
and if /i > 7a 2 , then the field points in the same direction for all |A| < a and the 
ground state is spin down. In this case, the reduced hamiltonian (18) is valid, and 
there exists a critical value of /i, fj,Q ~ 7a 2 

The net effect of the ip and x fields on the ground state is to shift the potential 
for A, and with proper scaling of 7, to shift the fermi energy. Let the potential 
V(X) be given by 

V(\) =(A 2 - X 4 /(2a 2 ))/a' -a<X<a 

(19) 

=00 |A| > a 

where a' is the Regge slope. This potential has maxima at A = ±a, and at 
these points V"(±a) = —1/a'. We have chosen the infinite wall to lie at the 
local maxima so that we can have a small magnetic field near these points and 
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still be able to compute the ground state. Adding the term in (18) shifts the 
maxima to ±a(l + 2^a' / ft) 1 / 2 and changes the second derivatives at these points 
to —(1 + 2 r ya ; / j3)/a' . Furthermore, if \i is tuned close to ixq then the magnetic field 
is small near the local maxima* 

The quartic potential is rather unwieldy to work with, so instead let us use the 
potential 

V(X) = - \ 2 /a < A < a 

(20) 

=oo A < or A > a. 

Now the effect of the boundary fields is to multiply V(A) by (1+270:'/ (3). The fermi 
energy for the potential in (20) is easily found using a semiclassical approximation 
for the density of states [2]. Inserting this approximation into an integral over 
phase space gives an expression for N, 



N = [3 J d\ J ^-6{e F -p 2 /2 + \ 2 /a'). 


The integral is straightforward, giving the relation 



(21) 



For our purposes we can approximate this equation as 

N = -^LU 2 + - A {-a'e ¥ ) log(-ep) + 0(e P )). 



(22) 



(23) 



27TVtt 

Shifting the potential leads to a new fermi energy e-p, which satisfies the equa- 



tion 

The value of \i does not effect this expression so long as the ground state remains 



N = (aVl + 27«V/5 + k-a/e F ) log(-e F ) + 0(e F )Y (24) 



* We should note that for 7 > 0, the points ±a(l + 2"/a'/0) 1 / 2 are in the region where the 
potential is infinite. But this won't matter because, as we will see later, aja' / (3 is much 
smaller than ((V(a) — e^/a') 1 / 2 . This is the distance from the wall where the potential 
equals the fermi energy. 
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spin down. Equating the expressions for N in (23) and (24) leads to the relation 

7 = ^2(eF-£F)log(-eF)- (25) 

It has been argued [4,5,19] that the proper double scaling limit is found by fixing 
(3y/a'ep to be constant while taking the limit N — > oo and (3/N to its critical value. 
Thus, we discover that 7 must have a logarithmic divergence in the scaling limit 
in order to shift (5ep by a finite amount. 

In order to understand the significance of this shift, let us consider the free 
energy in the scaling limit. Many authors have argued that this free energy is 
given by [3,4,5,19] 

-(pVa i e F ) 2 log(-£ F ) + \ log(-SF) 



(26) 



£(2^-1)^^)--). 



m=l 

where i?2m are the Bernoulli numbers. Based on this expansion, it seems clear that 
(— [3\fa'E~p)~ l should be interpreted as the closed string coupling, g cs . If we shift 
the fermi energy to ep and define k = (3(ep — Ep)/Va', then the relevant part of 
the free energy becomes 



47T/3 



— =( — ^-{(1 - a ! Kg cs ) 2 log(-e F ) - a? ng c 

(27) 



+ 2< QKfe > -% n(n -!)(„- 2)} 



1 1 

+ -log(-e F )--^(a / ^ cs ) n 

n—l 

_ V^/ 9 2m+l -i \ \B2m+2\ 2m (^ cs ) n (2m + n)\ \ 

£f >m{m + iy cs (2m)! )' 



Shifting the bulk cosmological constant shifts Ep and ep by an equal amount and 
thus keeps k fixed while scaling g ca . A surface of genus g with h holes should scale 
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as 



(fits) 2 *- 1 g ^ +h . Hence the terms in the expansion of the free energy have the 



proper scaling if k is interpreted as an open string coupling constant. 

If 7 < 0, then ep > ep and k > 0. In this case all terms in the expansion of the 
free energy are negative definite, except for those terms inside the curly brackets in 
(27). From (27) we also see that k = l/{g cs \fa i ) is a critical value for the coupling 
constant since ep — for this value, and thus the free energy for surfaces of genus 
g > 1 diverges when summing over the number of holes. The significance of this is 
not presently clear. 

Letting 7 > means that k < and therefore the terms in (27) alternate sign. 
Furthermore, /i must satisfy 



in order that the ground state is spin down. Therefore, in the scaling limit, the 
mass of the boundary generating fields diverges and hence, the dual graphs have 
Dirichlet boundary conditions. 

We can give a rough estimate of the correction to fiQ coming from the ferro- 
magnetic interaction. If we flip one spin up, then the ferromagnetic term will try 
to symmetrize the up spin among all the fermions. This is counteracted by the 
position dependent field which tries to push the spin to a point where the field 
is small. We will approximate these effects by assuming that the up spin is sym- 
metrized among the fermions less than a distance 2y from the top of the potential 
at A = 0. Hence, we will say that the average position of the up spin is at A = y 
and that it feels the 1/r 2 potential from the down spins that are greater than a 
distance y away from it. Thus, we estimate the shift in energy for turning up one 
spin as 



y 

where p(A) is the density of states, ^j. Plugging in the expression for 7 in (25) 
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(28) 




(29) 



and using the semiclassical estimate 



p(A) = ^L^\ 2 + e F a' (30) 
7rv a' 



gives 



Az? V2 Va*K\0g -£ F ) 2 2A /"0 

7r(3Va' pa z (3 

Minimizing AE leads to the estimates 

y ~ ? / i r vu/2' AjE ~ (^^-^^(^^(-eF)) - A/i ). (32) 

(a / /tlog(— 5f)J ' 

Therefore the correction to jiQ is of order 

A// ~ log(a'Klog(-e F )). (33) 

Actually, given the roughness of our estimate, it is not clear that we can distinguish 
between a double log and a constant. In particular, flipping the spin will certainly 
alter the density of states, and if the spin pushes the eigenvalues out far enough, 
then the double log should be replaced with a constant. In any case the correction 
Hq is much smaller than /iq 

Now suppose that 7 < and /j, 0. Let us calculate the energy to flip one spin 
at the bottom of the quartic potential. We can play the same game to estimate 
the critical value of /i. In this case turning up one spin shifts the energy by 

2y/2 f Vx 2 + a 2 2/i Kv / o 7 log(-£ F ) .2 



AE / dx ' +-hr~ % V ■ (34) 

PnvQJ x 2 (3 2a 2 K J 



y 

This then leads to the estimate 



fi ~ -(a/3)- 1 (-a'«log(-£ F )) 1 /3 (35) 

in order to keep AE small. Hence in order to have some relevance in the double 
scaling limit, the Neumann case also requires a divergent mass for the fundamental 
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fields, although the divergence is milder than in the Dirichlet case. In fact, the 
mass is small in the sense that //0/7 - ► in the double scaling limit. 

We would like to calculate the change in energy to turn up one spin if there 
are already s spins pointing up. This then leads to a density of up spins J^j. It 
is not clear how to do this, but a reasonable guess for this behavior is that the 
density scales as (3 log(— £p) f° r the Dirichlet case and (3 for the Neumann case. We 
postulate that the extra log term in the Dirichlet case is a result of the magnetic 
field being small at the top of the potential. Near this point, while there are fewer 
fermions around whose spin can flip, there are also fewer fermions to oppose such 
a flip. The guess is that the latter effect dominates the former. 

Turning to the one point functions of J n , it is clear that only (J™) is nontrivial. 
Actually, based on the structure of the hamiltonian in (17), the relevant operator is 
(1/A r )(tr0 n — 2J"), which couples to the fermion loops. Because of the structure 
of the vacuum, the one-point function for this operator is the same as the one- 
point function for (2/iV) tr n . Note that if \x » fio then those states that arc 
constructed by acting on the vacuum with J? decouple from the theory. In this 
case, the theory is essentially equivalent to a theory containing only closed strings, 
since now all J™ operators can be replaced by tr <p n operators. The only effect of 
tp and x i n this case is to renormalize the closed string coupling constant. 

One interesting operator is V = 1 — 2J®/N, where V generates linear transfor- 
mations of /i. If /i > fj,Q, then the one point function for V is 



In the case of dual graphs with Dirichlet boundary conditions, the length de- 



generates a shift of the boundary cosmo logical constant. For \x ^ fio, the boundary 
cosmological constant is, in some appropriate units, 



(V) = 2. 



(36) 



pendence of the free energy for a particular surface is (7//i) lensth . Therefore, V 



4(/i - no) 



(37) 



K 
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As in the case for closed strings, the cosmological constant should approach zero 
as the system nears the critical point. Since V generates transformations of order 
unity, the value of \x — /j,q should also be of order unity in order that it have some 
relevance in the scaling limit. In other words, if ji is shifted by order unity, then 
the energies of states with one spin pointing up are shifted by order 1//3, which 
is enough to have some effect on the physics. Thus, we find that /ib is of order 

(-log(-^F))- 1 . 

It is instructive to compare this with the behavior of the boundary cosmolog- 
ical constant for the one-matrix model. For the /c th multicritical point, the bare 
boundary cosmological constant scales as fi% = f3~ 2 /( 2k+1 ) /i R , where /iR is the 
renormalized value [12,13]. As k — > oo, c approaches 1 and the (3 dependence 
of /iB is (3 e with e — > 0. This result compares favorably with our result for the 
one-dimensional chain. 

The open string vertex operators should be constructed out of some combina- 

— * 

tion of the J n operators, which raise and lower the spins of the fermions. This 
suggests a possible connection to the Das-Jevicki collective field formulation of 
the tachyon field [6]. Recall that these authors showed there was a massless field 
associated with the fluctuations of the eigenvalue density, which they identified 
with the closed string tachyon. Analagously, it seems that the open string tachyon 
should be associated with the density of up pointing spins. This means that in the 
classical picture, there should be a connection between propagating tachyons and 
spin waves in the matrix model. 

In order to compute nontrivial open string scattering amplitudes, it is necessary 
to compute the spectra of states that have some spins up. It is not clear if it will 
be possible to do this exactly, but it is worthwhile to point out that the interaction 
term between the fermions in (17) is somewhat similar to the interaction studied by 
Calogero [20]. Calogero had the 1/r 2 potentials but not the spin-spin interaction. 
In any case, he showed that with a harmonic central potential, his theory was 
exactly solvable. It would thus seem not totally pointless to pursue exact solutions 
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of (17). Work on this issue and others is in progress. 
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